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The effective stress tensor of a homogeneous turbulent rotating fluid is anisotropic. This leads us 
to consider the most general axisymmetric four-rank "viscosity tensor" for a Newtonian fluid and 
the new terms in the turbulent effective force on large scales that arise from it, in addition to the 
microscopic viscous force. Some of these terms involve couplings to vorticity and others are angular 
momentum non conserving (in the rotating frame). Furthermore, we explore the constraints on the 
response function and the two-point velocity correlation due to axisymmetry. Finally, we compare 
our viscosity tensor with other four-rank tensors defined in current approaches to non-rotating 
anisotropic turbulence. 

PACS numbers: 47.27-i, 47.32.-y 

I. INTRODUCTION 

The properties and applications of rotating fluids constitute an important area of fluid mechanics [1]. In particular, 
the anisotropy consequent to the rotation has been a very important subject. For example, we have the classical 
Proudman- Taylor theorem, which says that, in the limit of fast rotation, the flow is so strongly anisotropic that 
actually becomes two-dimensional. Turbulence in the presence of uniform rotation, which is called rotating turbulence, 
is an example of anisotropic turbulence and an area of active research [1—4]. 

From a theoretical point of view regarding symmetry, the classical theory of fully developed turbulence assumes the 
maximal possible symmetry, namely, symmetry under translations and rotations, so it applies to ideal homogeneous 
and isotropic turbulence. However, in various situations such a high symmetry is not realistic and one must consider 
less symmetric turbulent states. The next most symmetric state is still homogeneous but the isotropy reduces to 
axisymmetry, that is, the full rotation group reduces to the group of rotations around a particular axis. The archetype 
of homogeneous turbulence with axisymmetry is rotating turbulence (naturally, the symmetry axis is the rotation axis). 

It was shown in Ref. [5] that perturbation theory for the randomly forced rotating Navier-Stokes equation generates 
anisotropic effective forces, in particular, the nondissipative force O x V 2 m. This suggests that one should find the 
complete set of allowed force terms. In this regard, it is useful to define the effective viscosity, which is a tensorial 
function of Jl (reproducing the known perturbative results for isotropic turbulence as f2 — > 0). Beyond perturbation 
theory (or any method of solution), our intention here is to study form first principles the consequences of axisymmetry 
in rotating turbulence. 

The possibility of anisotropy in the velocity correlation functions has been considered before in non-rotating fluids [6— 
12]. So, in these references, the anisotropy was attributed to other causes: existence of a mean flow or anisotropic 
forcing. In fact, in a homogeneous fluid the existence of mean flow effects, that is, the dependence of properties of 
the flow on its mean velocity, would contradict Galilean invariance. The flow can only depend on global kinematical 
features that involve accelerations, such as in a uniformly rotating fluid. A homogeneous but anisotropic forcing will 
induce anisotropy in the velocity field (the axial case is studied in Ref. [9]); but the physical origin of this anisotropic 
forcing and, therefore, the extent of the scale range affected by it are not clear. We think that rotating turbulence 
is a more natural example of anisotropic turbulence and with more physical applications. Moreover, this type of 
anisotropic turbulence has distinctive features (as was pointed out in Ref. [5]) because is an axial vector. Indeed, 
the force 17 x V 2 it or other terms of the same type would not be allowed if isotropy were broken by a polar vector, 
as in Ref. [9]. 

We remark that the characterization of the effective (or eddy) viscosity as a four-rank tensor has already appeared 
in the literature. For example, in Ref. [7] the authors show that a multi-scale method applied to the Navier-Stokes 
equation linearized with respect to a weak large-scale flow precisely produces an effective viscosity tensor if the basic 
fluctuating flow is not isotropic. However, as commented above, to determine the form of this tensor, one needs 
an explicit mechanism that breaks isotropy and preserves homogeneity. Otherwise, the basic assumptions and, in 
particular, axisymmetry, are not justified. 

On the other hand, since 17 is an axial vector, the effective viscosity tensor in rotating turbulence has distinctive 
features: for example, it has a pair-antisymmetric piece (which generates the above mentioned force) [5]. In addition, 
it will be shown here that a general treatment of the efective viscosity tensor in rotating turbulence requires new terms 
that couple to the vorticity or that are angular momentum non conserving (in the rotating frame) and, therefore, are 
forbidden in anisotropic non-rotating turbulence. 
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We shall first review the fluid equations in the rotating frame and the conditions for turbulence; we emphasize 
the transition from small scale isotropic turbulence to large scale anisotropic turbulence. Next, we introduce the 
viscosity in the standard manner [13] but without recourse to isotropy, which is replaced by only axisymmetry. All 
the components of the resulting four-rank tensor are determined with group theory arguments [14]. From this tensor 
we obtain the additional anisotropic force terms. Once seen the axisymmetry constraints on the viscosity tensor, 
we impose axisymmetry on the response function or (two-point) velocity correlations. In particular, the large-scale 
response function is related with the viscosity tensor. Finally, we try to connect the viscosity tensor with other 
four-rank tensors introduced in some current approaches to non-rotating anisotropic turbulence. 



II. EQUATIONS OF MOTION IN A ROTATING FRAME AND TURBULENCE 

The hydrodynamical equations for a fluid with density field p(x, t), velocity field u(x, t), pressure P(x, t) in a frame 
rotating with constant angular velocity ft are 

% + V-{pu) = 0, (1) 

+ (u ■ V)w = — VP - 2S1 x u - SI x (SI x x) + f , (2) 

at p 

where / accounts for an additional acceleration due to friction (which vanishes if diU = 0) and a homogeneous and 
isotropic external forcing, usually random (or periodic as in Ref. [7]), which serves for keeping the total kinetic energy 
constant. 

We assume that the fluid is incompressible, with constant density, so the continuity equation becomes V • u = 0. 
So if we define p = P/p every reference to the density disappears, and we have two equations for the two unknowns 
u and p. To solve for p, it is useful to separate Eq. (2) into independent longitudinal and transverse equations. Since 
u = ul + ut (such that V x ul = V • ut = 0) and ul identically vanishes, the longitudinal equation becomes just 
a constraint relating p with spatial derivatives of u, namely, 

p=^(fixi) 2 -^ [di (ujdjUi) + 2eijkSljdiUk] ■ (3) 



Solving for p, the equation for u = ut is 



Oil 

— + V[(u • V)u] = -V(2Sl x u) + f , (4) 

where the projection operator V onto transverse (or solenoidal) fields is given by 

P = 1-V^V-. (5) 

In Eq. (4), if u is tranverse so is / and vice versa. We call Eq. (4) the transverse rotating fluid equation. If we 
substitute for / an isotropic viscous force, it becomes the transverse rotating Navier-Stokes equation. Note that the 
transverse rotating fluid equation (4) is translation invariant (assuming that / is homogeneous), in contrast with 
Eq. (2). Therefore, its solutions are homogeneous velocity fields and, furthermore, one can make use of the Fourier 
transform. 



A. Homogeneous rotating turbulence 

The homogeneous rotating turbulent state is defined by a velocity field with large fluctuations but such that the 
mean velocity is negligible in the rotating frame. Let us see how to characterize this state in terms of nondimensional 
parameters and how it is related with the homogeneous and isotropic turbulent state. 

Since we have the rotation velocity as additional parameter, we can define two nondimensional parameters, namely, 
the Reynolds and Rossby numbers. While the Reynolds number Ro — UL/v measures the relative importance of the 
nonlinear and viscosity terms in the Navier-Stokes equation, the Rossby number Ro — Uj (LSI) measures the relative 
importance of the nonlinear and Coriolis terms in the rotating Navier-Stokes equation (U is a reference velocity or 
the variation of the velocity over the length L that characterizes the system size). In principle, Ro ^> 1 indicates that 
rotation effects are negligible and, viceversa, Ro <C 1 indicates that they are dominant. However, the latter condition, 
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namely, dominance of rotation effects over nonlinear and viscous (and even dynamic) effects leads to the situation in 
which only the linear Coriolis force is relevant, giving rise to extreme two-dimcnsionalization of the flow (as in the 
Proudman- Taylor theorem) but without turbulence, ft is necessary that the two numbers Ro and Re play a role in 
specifying the regime of interest, that is to say, the regime with rotation effects (Ro » 1) and turbulence (Re » 1). 
Or one may introduce the Ekman number Ek = v/ (HL 2 ) (in addition to Ro), which is the ratio of the Rossby number 
to the Reynolds number and measures the relative importance of the viscosity and Coriolis terms [1]. Then one must 
demand Ek <C Ro in addition to Ro <C 1. 

To clarify the preceding condition, let us consider relevant length scales. First, let us recall the role of the dissipation 
scale. In ordinary homogenous and isotropic turbulence, K41 theory [8] makes the dissipation rate per unit mass e 
the basic quantity and introduces the dissipation scale A = (i^ 3 /e) 1//4 . Using fl instead of v, we can form with e the 
length scale I = (e/S! 3 ) 1 / 2 . If we begin with small £1 (for fixed e) such that I » L, rotation effects must be negligible 
all over the fluid system of characteristic length L. Therefore, the precise condition for neglecting £7 is £1 <C (e/L 2 ) 1 / 3 
(equivalent to I > L). Given that (£/A) 4 / 3 = Re > 1, the parameter (e/L 2 ) 1/3 /^ = v / (£IL 2 )(L / = EkRe = 
Ro, so the condition for neglecting f2 is just &> 1. As £7 grows and, therefore, I diminishes such that I < L, 
rotation effects become appreciable. We then have one scale range with rotating turbulence, namely, between i and 
L, and an othe r with isotropic turbulence, namely, between A and I. The latter range holds as long as A < I, that 
is, £1 < \J ejv. As £ becomes smaller than A, the rotation effects dominate over the nonlinear effects and the flow 
becomes strongly two-dimensional. 

The interesting values of £1 are such that there are the two scale ranges, respectively, with isotropic turbulence on 
small scales and anisotropic turbulence on larger scales. Of course, this happens when A « ^ < £. Then the viscosity 
or correlation functions on scales between A and I are essentially isotropic whereas the effective viscosity or correlation 
functions on scales between I and L are axisymmetric. 



III. THE AXISYMMETRIC EFFECTIVE VISCOSITY TENSOR 



To introduce the viscosity tensor, it is convenient to follow the general reasoning [13] which starts by writing the 
fluid equation in local conservative form, as 

dui dXVij 

P ~dt = ~dx~' ij = ~ pUlU 3 + Ti > ' ( 6 ) 
and finds the deviatoric part of the stress tensor 

Ti^-PSii+oa (7) 

due to internal relative motion (viscosity) from general principles. The first principle is that the velocity gradient is 
small, which allows one to consider only first derivatives of the velocity. Next, the viscous stress tensor ct^ is taken 
proportional to the velocity gradient and, furthermore, its antisymmetric components (vorticity) are excluded, so 
that the stress is proportional to the rate of strain u mn — dr m u n ^ — (du n /dx m + du m /dx n )/2 (this characterizes 
Newtonian fluids) [16]. The following crucial assumption is isotropy, which leads to the existence of only two propor- 
tionality constants (shear and bulk viscosities). As we cannot make this assumption here, we are left with just the 
proportionality relation 

Vijmn ^mn ? (8) 

like in the analogous relation in the theory of elasticity that expresses that the stress is proportional to the strain [15]. 
Therefore, the symmetry properties of the tensor T]ij mn , which we call the "viscosity tensor", are similar to the ones 
of the elastic modulus tensor, namely, symmetry under exchange of indices within the first and second pairs of indices 
and, in addition, symmetry under exchange of the first and second pairs of indices (pair symmetry). However, we shall 
further allow for pair antisymmetry; namely, we write Vijmn as a sum of a pair-symmetric (S) and a pair-antisymmetric 
(A) part [5]: 

Vijmn — -^(l]ijmn + Vmnij) + ~^(Vijrnn ~ Vmnij) — Vijmn + Vijmn • (9) 

So, generically, the "viscosity tensor" has 36 independent components, of which 21 belong to the pair-symmetric part 
Vijmn an d 15 belong to the pair- antisymmetric part Vijmn- 

The axial symmetry of the equations of motion reduces the number of independent components of both vfjmn 
and Vijmn- The 21 components of the generic pair-symmetric tensor can be divided into two sets with 15 and 6 
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components, respectively, the former corresponding to the totally symmetric tensor. The respective components are 
constructed in the appendix as linear representations and called 15S and gS. Further imposing axisymmetry, the 
pair-symmetric tensor can be constructed from fij and Sij as 

Vijmn — ^li^ij^mn ~t~ SimSjn ~t~ Si n Sj m ) 

~\~ ^2(^i^j^mji ~i~ Qm^>nSij ~t~ Qi^m&jn ~i~ ^j^m^iro ~t~ ^i^n^m ~i~ ^j^n^im) 
-|- d^Qi^lj^lm^lji 

+ CliSijSmn 

+ a 5 (ttittj5mn + Qm^lnSij) . (10) 

There are five independent components, to which we attach scalars oi, . . . , ag (which can depend on ft 2 ). In comparison 
with the form given in Ref. [5], this expression has been arranged so that the three first tensors (with coefficients 
a ii a 2i a 3) are totally symmetric in their indices. 

The generic pair-antisymmetric tensor has 15 components, constructed in the appendix as the linear representation 
15S". The pair-antisymmetric tensor with axisymmetry needs, in addition to and Sij, the totally antisymmetric 
tensor e^-fc and is 

Vijmn — ^l^qi^qimSjn ~t~ ^qin&jm ~\~ ^qjm&in ~t~ ^qjn&im) 

-\- b2^lq^qim^j^n ~\~ ^qin^j^m ~f" ^qjm^i^n ^qjn^i^m) 

+ b3(SliSlj5 mn - n m £l n 8ij) ■ (11) 

We observe that the axisymmmetry has reduced the number of independent components from 21 to 5 for the pair- 
symmetric part and from 15 to 3 for the pair-antisymmetric part. This reduction can be explained by considering the 
reduction of linear tensor representations under rotations (see appendix and Ref. [10]). The reduction under rotations 
is performed by extracting traces, which are rotation invariant but not linear invariant. There exists a canonical 
procedure for doing this trace extraction [14] but we can clarify the procedure by noting that the properties of the 
expressions in Eq. (10) or Eq. (11) under rotations are determined only by the vector Hi (Sij and e^fc are rotation 
invariant). For example, the terms with coefficients a\ and 04 clearly correspond to scalars (J = 0), the only terms 
allowed by isotropy. Furthermore, an expression with the tensor product of n O's corresponds to the representation 
J = n, usually, with an admixture of lower J representations. Therefore, each coefficient corresponds to a definite J 
representation, but, in order to obtain the correct tensorial expression of each representation, we need to remove the 
lower J representations by extracting traces. This induces a linear redefinition of the coefficients within each linear 
representation. 

Finally, we remark that the terms in Eq. (11) with coefficients b\ and 6 2 would not be allowed if isotropy were broken 
by a polar vector, because the respective terms would be odd under parity (in general, the parity of the representation 
J associated to a polar vector is (— ) J ). 



A. Traceless components and incompressibility constraint 

The preceding tensorial expressions for the viscosity have a part that couples to the velocity divergence uu. More- 
over, they give rise to an isotropic part of the viscous stress tensor Uij (that is, proportional to Sij). Therefore, the 
viscosity tensors must be further decomposed into traceless and trace parts. For incompressible flow we only need the 
traceless components such that r\ijkk — Vkkmn = 0. They can be extracted by subtracting traces from either Vijmn 
or vtjmn- I n the general case, that is, with no axial (or any other) symmetry, those tracelessness conditions remove 
6 + 6—1 = 11 components (the condition r/kkii — appears twice), leaving 25 components. To be more precise, the 
conditions vfjkk = remove 6 components of vfjmn (with J = 2, 0, corresponding to 6 5), and the conditions rjf- kk = 
remove the 5 components of rj^ mn corresponding to J = 2. 

Indeed, a straightforward calculation yields: 

Vijmn ~ T^VijkkVmn ~ T^Vkkmn^ij ~gVkkll^ij^rnn = a l\VimVjn ~\~ OinOjm ~ T^VijVmn) H~ 

4 4 

&2 [^i^m^jn ^j^m^in H~ ^i^n^jm ^j^n^im "(Hjflj^mn ^m^n^ij ) H~ ^ij^mnj 
Vijmn ~ ^Jlijkk^rnn ~ T^Vkkmn^ij ~ ^l^q( € qim^jn H~ e qin&jm H~ ^qjrn^in H~ £qjn$im) 
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~t~ ^2^g y^qim^j^n ~i~ ^qin^j^m ^qjm^i^n ~i~ ^gj'n^z^m)* (^^) 

The number of coefficients has been reduced to three for r]fj mn , corresponding to the J = 4, 2, representations, and to 
two for r]{j mn , corresponding to J = 3, 1. It is natural that they together constitute the Clebsh-Gordan decomposition 
of the tensor product of two J = 2 representations (with dimension 5x5 = 25) [14]. 

There is another set of tracelessness conditions, namely, -q^mj = 0, but there is no physical reason to impose them. 
However, note that the six conditions ~qfj m j = remove the J — 2, representations, just leaving J = 4, while the 
three conditions vjjmj = remove the J — 1 representation, just leaving J = 3. Therefore, this last set of tracelessness 
conditions would select the highest J representations, corresponding to the coefficients 03 and 62- 



B. Viscosity tensors with antisymmetric pairs 



Two crucial assumptions in the reasoning at the beginnning of Sect. Ill are that the viscous stress tensor is symmetric 
and that it does not depend on the vorticity (the vorticity tensor is = duUj] = (diUj — djUi)/2). They lead to 
a viscosity tensor with symmetry under exchange of indices within the first and second pairs of indices (symmetry 
by pairs). Those two assumptions are commonly accepted since they are based on basic physical principles: on the 
one hand, the stress tensor can always be chosen symmetric because of angular momentum conservation; on the 
other hand, a uniform rotation (leading to a constant vorticity) cannot induce stresses, so a dependence of the stress 
tensor on vorticity is forbidden. However, both principles, namely, angular momentum conservation and absence of 
stresses in uniformly rotating fluid, fail in a rotating frame. Therefore, we are allowed to consider viscosity tensors 
with antisymmetric pairs of indices. We have three types: (i) tensors Xijmn with the first pair symmetric and the 
second antisymmetric, which account for an angular momentum conserving coupling to vorticity, (ii) the symmetric 
type £,ij mn , that is, tensors with the first pair antisymmetric and the second symmetric, which account for an angular 
momentum non-conserving coupling to strain rate, and (iii) tensors dj mn with both pairs antisymmetric, which 
account for an angular momentum non-conserving coupling to vorticity. 

The most general tensor with the first pair symmetric and the second antisymmetric has 18 components (see the 
appendix). Its axisymmetric form is 

Xijmn {ci^ij "i" C2^i^j j^imu^i ~t~ 
C-Z^S^i^m^ jn ~t~ ^j^m^in ^i^n^jm ^j^n^tm) ~^ ^A^imn^j ~t~ ^jrarS^i) • (1^) 

The constants ci,C3,c 2 correspond to J = 1,2,3, respectively, forming the linear representation i5 SA, whereas C4 
corresponds to J = 1 and 3SA. Imposing that the tensor be traceless in its first two indices, that is, Xamn = 0, 
relates the coefficients c\ and c 2 (the tensor is automatically traceless in the second pair of indices). Therefore, the 
traceless tensor contains the J = 1,2,3 representations, corresponding to the Clebsh-Gordan decomposition of the 
tensor product of the 3 = 1 and J = 1 representations. 

There is an analogous axisymmetric structure for the symmetric type £r/ mrl , involving 15 AS and 3 AS, and with 
coefficients d-y , . . . , C4. 

Finally, the tensor Qj mn with both pairs antisymmetric has 9 components, which the axisymmetry reduces to 

Cijmn — ^l(^im^ jn $in& jni) ~l~ 
d^i^i^m^jn Qj^m^in ^i^n^jm ^j^n^im) ~t~ ^S^imn^j ^jmn^i) • (1^) 

The constants <ii,G?2 correspond to J = 0,2, respectively, forming the representation qA (which is pair symmetric), 
whereas (I3 corresponds to J = 1 and 3A (which is pair antisymmetric: even though it may not seem obvious, 
timn^j — tjmn^i = —^mij^n + ^nij^m)- The tensor defined by Eq. (15) is trivially traceless in both pairs of indices 
and corresponds to the Clebsh-Gordan decomposition of the tensor product of two J = 1 representations. 

C. Effective forces associated with the viscosity tensor 

The total viscosity tensor T = rj + x + £, + (is defined by 

&ij ^ijmn ^m^n • (16) 

The force derived from this stress tensor is 



fi d&ij / dxj — Tijmn ^jm^n • (l^) 
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The expression that results by substituting the full axisymmetric expression of ry TO „ is fairly complicated: suppressing 
gradient terms, we obtain 

/ = (01 - rfi)v 2 w - 61 (n x v 2 «) - 6 2 (n • v) 2 (o x «) - (6 2 + 4) (n • v)(o x v)(n • «) 
+ (6 2 + c 2 ) n (n • v)(n-w) + (C4 + 4 + d 3 )(n • v)w + (a 2 - c 3 + 4 - d 2 )fi v 2 (o • u) 

+ (a 2 + c 3 -c 3 -d 2 )(rJ-V) 2 M + a 3 J7(rj- V) 2 (fl-M). (18) 

Several remarks are in order. Note that the fifth and sixth terms of the force involve the vorticity u> = V x u and are 
proportional to an odd power of fi. The terms preceding them are also proportional to an odd power of fi, except the 
first one, which is isotropic. The remainig three anisotropic terms, which neither involve the vorticity nor any vector 
product, are equivalent to the anisotropic force written in Ref. [9]. If we had considered only the tensor r] S to derive 
the force, we would have obtained precisely these three terms but the first couple of them would have had the same 
coefficient (a 2 ). As we use the complete tensor r we have instead that some coefficients are redundant: inspecting Eq. 
(18), we see that there are two redundant coefficients among C4, c 4 , d 3 , two redundant coefficients among a 2 , c 3 , c 3 , d 2 , 
and one redundant coefficient among a\,d\. 

After taking into account that V-u = and suppressing gradient terms, only remain the coefficients of the part of r 
that is traceless in the first and second pair of indices. Gradient terms are longitudinal and the physical force must be 
transverse (solenoidal); but, after removing these terms, the force is still non transverse and must be projected with 
the nonlocal operator V of Eq. (5). This operation brings back two supressed gradient terms, namely, V(fi • oj) = 
and V[(f2 • V)(fi • u)], in addition to producing nonlocal gradient terms. 

Finally, we remark that all the terms in Eq. (18) coming from odd- J components of r, that is, the ones with odd 
powers of fi (with coefficients 61, 6 2 , c 2 , c' 2 , C4, c 4 and d 3 ), would not be allowed if isotropy were broken by a polar 
vector. 

D. Dissipation 

The dissipated power is 

— J d 3 x u • f = — J d 3 x Ui djOij = j d 3 x djUi o~ij = 

J d X UijTjij mn U mn -\- J d X UijXijmn^rnn J d X ^ij^ijmnUmn J d X tOijC,ijran^ran , iX^) 

where we have assumed that the velocity vanishes on the boundary to remove the surface integrals, that is, we have 
assumed that there is no work made by external sources. 

As remarked in Ref. [5], r\ A does not lead to dissipation; neither does ( A . Moreover, if Xijmn = £mnij, the respective 
terms cancel in Eq. (19). All these nondissipative components of r do not properly belong to the viscosity tensor, 
although they give rise to forces with dynamical effect. On the other hand, since the dissipation cannot be negative, 
we can deduce some positivity conditions on the proper coefficients of the viscosity tensor: a\ > 0, —d\ > 0, etc. 



IV. AXISYMMETRIC FORM OF THE RESPONSE FUNCTION AND VELOCITY CORRELATIONS 



It is useful to study the symmetry constraints on the response function and velocity correlations. Here wc determine 
the most general axisymmetric forms of these quantities in the small-wavenumber limit (corresponding to large-scale 
features). The theory of axisymmetric tensors has been the subject of previous analyses of anisotropic turbulence; 
in particular, it has been treated in papers by Chandrasckhar [6] and by Arad et al [10]. The former uses the 
old formalism of invariant theory whereas the latter uses the theory of group representations. Unfortunately, both 
consider only the application to correlation functions in real space, while we are interested here in correlation functions 
in Fourier space (spectral functions). Therefore, the theory of axisymmetric tensors as is developed in those references 
must be adapted to Fourier space. Actually, the spectral two-point velocity correlation function in rotating turbulence 
has already been studied by Cambon and Jacquin [2] and we shall use their results. 



A. Axisymmetric form of rank-two tensors 

We consider a second-rank tensor that depends on the wave vector k (since we use Fourier space), in addition to 
the angular velocity ft. The general form of such a tensor as a linear combination of the tensorial products k{kj, 
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fijfij, fc,f2j, rijfcj and the unit tensor <5jj is 

T i:j (fc) = Afc.fcj + BOiOj + CTcjOj + + , (20) 

where A,B,C,C',E are arbitrary functions of k, fl and fe • fl. However, a more general expression results upon 
introducing the unit antisymmetric tensor or, equivalently, the vector n = k x fl (assuming that k and fl are not 
parallel) and the corresponding tensor products: 

Tijiyk^j = Ak^kj ~\~ Cki^lj ~\~ DkiTtj 

+ D'mkj + F'mflj + Gn.rij . (21) 

This expression with nine coefficients is the most general one, because any vector (to be included in a tensor product) 
can be expressed as a linear combination of fe, fl and n. 

We remark that expression (20) corresponds to the ordinary quadratic form of Ref . [6] , where the terms with the 
unit antisymmetric tensor are named "skew" forms. This name refers to its reflection (or parity) character: if the two 
vectors employed in the tensor products of the ordinary quadratic form are polar, this form is parity invariant (even 
parity), whereas the skew forms change sign under reflections (odd parity) since the vector product is axial. In our 
case, we begin with a polar vector k and an axial vector ft, so their vector product n is polar. Hence, the terms of 
Eq. (21) that change sign under reflections are the ones with only one ft. 

Instead of the basis formed by fc, fl and n, it may be more convenient to use an orthonormal basis. Any couple 
of linearly independent vectors determine an orthonormal basis, in particular, the two vectors fc and ft lead to the 
one given by fe/fc, = fc x ft/ \k x ft\ and = k x / |fc x eW| [2]. Note that the vector ft is axial, so 
is polar but is axial. Interchanging the role of fl and fc, wc get a different orthonormal basis, with the vector 
in common: e 3 = ft/fl, e 1 = fc x ft/ |fc x ft\ = and e 2 = fl x e 1 / \ ft x e 1 ]. This basis (which we denote by 
superindices without parentheses) is more adequate. Note that both e 1 and e 2 are polar. 

Any rank-two tensor can be expressed in the latter basis as 

Ta = f pq e\e) . (22) 

There are three pieces in that are independent under rotations: the trace, the antisymmetric part and the 
traceless symmetric part. This is the Clebsch-Gordan decomposition of the vector tensor product into the irreducible 
representations J = 0, J = 1 and J = 2 of the rotation group. However, to classify the behaviour of the components 
of the second-rank tensor under rotations around fl, that is, under the two-dimensional rotation subgroup 0(2) of 
the full rotation group 0(3), it is best to use the given basis (components T pq ). The irreducible one-dimensional 
representations of 0(2) are complex, labeled by an integer M {—J < M < J). The real irreducible representations 
are labeled by \M\ and are two-dimensional (except the scalar M = representation) [14]. We have that e 3 is a 
scalar, and {e^e 2 } form the real \M\ = 1 representation. Consequently, T33 is the scalar M = representation, 
^13, T23, T31, T32 belong to two \M\ = 1 representations, and the remaining components in the 2x2 block matrix can 
be subdivided into its trace (M = 0), its antisymmetric part (M = 0), and its faceless symmetric part (\M\ = 2). 
Furthermore, it is not difficult to ascribe each M representation to a definite J representation. 

B. Axisymmetric form of the response function 

The response function is defined by 



(23) 



5fj(k,uj) 

(introducing a non-random part in the external forcing /). So we can write, at linear order in /, 

(u i (k,u))=G ij (k,u)f j (k,u) . (24) 

Conversely, 

GT j 1 (k,uj)(u j (k,ij))=f i (k,uj), (25) 
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which tells us, on account of Eq. (17), that the quadratic term in the expansion of G ij 1 (fc, 0) in powers of k is related 
with the viscosity tensor. To be precise, we have that 



1 d Gjj 
g ijmn ■- 2 dkmdkn 



2 (Timjn ~t~ Tinjm) i ( 2 ^) 

fe=0 



defining what we call the (four-rank) response tensor. This tensor has symmetry in the pair inn, so it has 54 
independent components, while t (in the generic case) has 81 components. Indeed, the 27 components of the tensor 
\{jimjn — Tinjm) do not contribute to the response function. 



1. Mapping the viscosity tensor to the response tensor 



We can take in Eq. (26) 77, x, £ or Q for r. On the other hand, gij mn can be decomposed into zj-symmetric and 
ij-antisymmetric parts, corresponding to the respective parts of the response matrix. Therefore, gijmn has 6 x 6 = 36 
components with symmetry in both pairs (belonging to the £ representation) and 3 x 6 = 18 components with 
antisymmetry in the first pair and symmetry in the second pair (belonging to the AS representation). 

Let us first analyze the components of g coming from rj s . We note that gij mn = (vfmjn + Vinjm)/ 2 = {Vjnim + 
^fmm)/ 2 = gjimn- Furthermore, this tensor is pair-symmetric: g mnij = {rj minj +V m jni)/ 2 = (vLjn +vf mm )/ 2 = 
(vfmjn ^Vinjm)/ 2 ~ gijmn- So the 21 pair-symmetric components of n s (representations 15S and 65*) are transformed 
by Eq. (26) into the 21 pair-symmetric components of g; in particular, the totally symmetric representation i 5 5 
of i] s is left invariant. Given that we can substitute i] S by ( s in the preceding equations, we conclude that the 6 
pair-symmetric components of £ s (qA) are transformed by Eq. (26) into 6 pair-symmetric components of g (linear 
combinations of 15 S and qS). 

We also note that if we take x f° r T m Eq. (26) and symmetrize in ij, the tensor gij mn — \{ximjn + Xinjm + 

Xjniin 1 Xjnim 

) is pair antisymmetric, owing to the symmetry of x- An analogous property is fulfilled by the tensor 
gijmn constructed in the same way from £. So the 15 components of \ from 15SA or the 15 components of £ from 
15AS are transformed by Eq. (26) into the 15 pair-antisymmetric components of g (representation 15S"). On the other 
hand, it can be proved that the x or ^ belonging to representations 3SA or 3^5*, respectively, yield vanishing g (they 
contribute instead to T im j n ~ T injm)- 

As regards the ij-antisymmetric part of g ijmn , note that g ijmn = {vtm 3 n+ r ltnjm) / 2 = (^Im-ljim)/ 2 = -ffjimn- 
So the 15 components of -q A (representation 15S") are transformed by Eq. (26) into 15 components of gtjmn with 
antisymmetry in ij (15 AS). Given that we can substitute r\ A by ( A in the preceding equations, the remaining 3 
components of ( A (3A) are transformed into three components of gijmn forming the representation 3AS. Finally, 
we also have the mapping SA -> AS given by g ljm n = jiXimjn + Xin 3 m - Xjmin - Xjnim) and a similar mapping 
AS -» AS. 



2. Axisymmetric form of the response tensor 



The preceding mapping has been established with full generality, without considering any particular spatial sym- 
metry. If we take axisymmctry into account, Eq. (26) provides the linear relations between the coefficients in the 
axisymmetric form of gijmn and the coefficients in the axisymmetric forms of n, x, £ or C The 54 components of 
gijmn belong to the S and AS representations, therefore, gijmn has 8 + 4 = 12 coefficients: the axisymmetric form of 
the components gijmn with symmetry in the first pair is like the forms of rj in (10) and (11), with other coefficients, 
say ct\, . . . , 0:5 , /3 2 , /?3; the axisymmetric form of the components gij mn with antisymmetry in the first pair is like 
the form of £, with other coefficients, say 71, ... , 74. 

The above-mentioned coefficients in the response tensor can also be obtained by expanding in powers of k the 
axisymmetric expression of Gr. 1 (21). Considering that the unit antisymmetric tensor does not appear in the tensors 

with coefficients a\, . . . , as, 03, 73, the corresponding part of G" 1 is given just by expression (20). Then the coefficients 
ai, . . . , as, /?3, 73 must arise by expanding A, B, C,C',E in powers of k such that the total expression is of second 
degree in k. In particular, B = B\k 2 + B 2 {k ■ fi) 2 and E = E\k 2 + E 2 {k ■ 1~2) 2 , while A is a constant and the 
coefficient functions of the terms that are of first degree in k, namely, C,C, can only be expanded up to the first 
order (proportional to k-Sl). Therefore, this expansion just doubles the coefficients of Qjfij and 5{j, producing seven 
coefficients altogether. We can divide them between the six coefficients arising from the symmetric (C = C) 

and the one corresponding to the antisymmetric GfT 1 ,, namely, C — C . 
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The part of that includes the unit antisymmetric tensor, with coefficients D,D',F,F',G in Eq. (21), corre- 
sponds to /3i , fa , 7i > 72 and 74. They can be divided into /3i,/32 for the symmetric and 7i>72 and 74 for the 
antisymmetric G^.j. 

3. Higher-rank axisymmetric response tensors 

An expansion of G,^ 1 in powers of k to orders higher than the quadratic order yields response tensors similar to 

9ijmn in Eq. (26) but of higher rank. Given that G" 1 must be parity symmetric, only even powers of k can appear. 
For example, the next higher-rank response tensor gij mnp q is symmetric in the last four indices and, therefore, has 
9 x 15 = 135 components, but this number is reduced by the axisymmctry. The following higher-rank response tensors 
are progressively more complex, of course. 



C. Axisymmetric form of the two-point velocity correlation 

Let us introduce the spectral two-point velocity correlation: 

{u i {k,oj)u j {k',J)) = (2ir) 4 U t3 (uj,k)5(uj + uj')5 3 (k + k'). (27) 

We have that Uij(cJ,k) = Uji(—w,—k). Furthermore, transversality implies that kiUij — hUji = 0. So, in the 
isotropic case, the spectral two-point velocity correlation is given in terms of only one function: 

Uij(u, k) = V tl {k)U{uj 1 k). (28) 

Taking into account that equal-time correlations are more useful, let us define 

UiAk) = j '^Uiji^k) (29) 

(assuming that the integral is convergent), so that 

(Ui(fc,t)Mj(fc',t)) = (27r) 3 ^(fc)5 3 (fe + fc'). (30) 

As demonstrated in section IV A, the general axisymmetric rank-two tensor has nine independent coefficient functions, 
but the transversality conditions reduce their number. The number of independent conditions is five, so just four 
coefficient functions remain independent, namely, the ones corresponding to the tensor products of the transverse 
vectors e^ 1 ) and e^ 2 \ It is convenient to use the basis corresponding to circular polarizations TV = — ie^ 2 \ 
N* = + ie^ 2 \ so that the resulting tensor can be written as [2]: 

Uij(k) = e(k)V tJ + XlzWNiNj] + ih(k)eiji^. (31) 

The quantities e(fc) and h(k) are the energy and helicity spectrum and z(k) is a "complex deviator". They all 
are even functions of k. The preceding form is equivalent to the form with the products ef^e^p, on account that 

NiNj = e«e« - ef >ef - i(e«ef + ef ef), V l3 = ef^ef + ef ef , and e ijlkl /k = ef^ef - efef. 

Velocity correlations for more than two points lend themselves to be expressed in similar though more complicated 
ways. 



V. CONNECTION WITH SOME APPROACHES TO ANISOTROPIC TURBULENCE 

We have already mentioned that fourth-rank tensors associated to anisotropic turbulence (but without rotation) 
have been studied before; for example, in Ref. [7]. More recently, in Ref. [11] has been defined one fourth-rank 
tensor for a flow with a constant strain rate. Another fourth-rank tensor is defined in Ref. [12] in connection with 
the linearization of a closure equation in the presence of weak anisotropy. We now explore connections between our 
anisotropic viscosity tensor and those tensors. 



10 



The fourth-rank tensor Cij mn (k) of Ref. [11] expresses proportionality between the contribution to the correlation 
Uij(k) (defined in Eq. (29)) from anisotropy and the constant strain rate producing the anisotropy: 

&lAij{l%) = Cijmn{k) ^>mn > (32) 

where the strain rate u mn is constant. The Reynolds stress tensor 

/d 3 k 
— U ij {k)=U ij (33) 

has a deviatoric part that, according to Eq. (32), is proportional to the strain rate, the proportionality constant being 
the integral of the tensor Cij mn (k): 

d3k RU (V\ [ d3ft 



Mij= / 7KZX3 Su n( k ) = / C ijmn (k) u mn . (34) 



It is also possible to assume that the Reynolds stress tensor (33) and the strain rate have some mild dependence on 
the spatial coordinate x. The corresponding generalization of Eq. (34) is a phenomenological (mean- field) closure 
relation that can be justified with a multi-scale method applied to the Navier-Stokes equation linearized with respect 
to the a;-depcndent (large-scale) mean flow [7]. Comparing this mean- field relation with Eq. (8), we deduce a relation 
between our viscosity tensor rjij mn and Cij mn (k), namely, 

Vijmn — P j ^2^3 Cijmn(k)- (35) 

We must note, however, that the form of Cij mn (k) in terms of projectors Vij(k) proposed in Ref. [11] leads to the usual 
isotropic rjij mn . Indeed, one needs an additional quantity such as the vector f2, to define an anisotropic viscosity. 

More sophisticated closure schemes involve relations between the three and two-point velocity correlation functions. 
The Navier-Stokes equation leads to an equation involving these two types of correlations, first derived by von 
Karman and Howarth assuming isotropy. Chandrasekhar [6] developed a theory of axisymmetric tensors to generalize 
this equation to axisymmetric turbulence. As remarked by Frisch [8], it is easy to derive a fully anisotropic version 
of the Karman-Howarth equation, which he calls the Karman-Howarth-Monin equation. In Ref. [12], the Fourier 
transform of this equation is used as the basis of a closure scheme related with the direct interaction approximation 
(DIA), in which the function Uij{k) satisfies (in stationary conditions) a nonlinear equation: 

D m n(k) = I mn (k) - vk 2 U mn (k), (36) 

where I mn is an integral operator quadratic in Uij . In addition, we have introduced an external random forcing, absent 
in Ref. [12], which is Gaussian and white in time and, hence, is represented by the spectral two-point correlation 
Dmn(k) (see Ref. [8] for a general description of closure equations). If the molecular viscosity v vanishes, the external 
forcing is not necessary as in Ref. [12]. However, the introduction of D mn allows us to substitute the four-rank tensor 
defined in Ref. [12] by a four-rank tensor more useful to connect with the viscosity tensor. 

If the forcing is isotropic, we expect that Eq. (36) has isotropic solutions. One can then linearize this equation 
around an isotropic solution, namely, 



oU mn {q) = 



5U t] {k) 



5U i3 {k), (37) 



where SD mn represents an anisotropic perturbation of an isotropic forcing such that D mn — DV mn (note that isotropy 
implies that Uij = eVij, according to Eq. (31)). The solution of this linear equation is obtained by inverting the 
matrix of pairs of indices, deriving a sort of tensorial response function, 



G (ka)- 5Ulj{k) 



6D mn (q) 



(38) 



which measures the response to the anisotropic perturbation. Considering the role of the molecular kinematic viscosity 
v in Eq. (36), we can tentatively define a kinematic viscosity tensor as 



(39) 

This relation between a tensorial viscosity and a response tensor is an alternative to Eq. (26), valid when we replace 
the original Navier-Stokes equation with the closure Eq. (36). However, the actual computation of Vij mn necessarily 
leads to an isotropic tensor, since there is nothing in Eq. (38) capable of breaking rotation invariance. 
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VI. CONCLUSIONS 

We have applied symmetry principles to homogeneous turbulence subjected to uniform rotation, focusing on the 
four-rank tensor defining the linear relation between the stress tensor and the velocity derivatives, which we call the 
"viscosity tensor". The most general tensor comprises five parts: 

• a tensor -q symmetric by pairs of indices and pair symmetric, accounting for the usual proportionality relation 
between the (anisotropic) stress and the strain rate; 

• a tensor r/ A symmetric by pairs and pair antisymmetric embodying a new relation between the stress and the 
strain rate, typical of rotating fluids, since it does not lead to dissipation; 

• a tensor \ symmetric in the first pair of indices and antisymmetric in the second, which accounts for a stress 
tensor coupling to vorticity; 

• a tensor £ antisymmetric in the first pair of indices and symmetric in the second, which accounts for the 
antisymmetric part of the stress tensor (angular momentum non-conserving) that couples to the strain rate. 

• a tensor ( antisymmetric in both pairs of indices which accounts for the antisymmetric part of the stress tensor 
that couples to the vorticity. This tensor can be further decomposed into pair-symmetric and pair-antisymmetric 
parts, like 77. 

Group theory helps us to find the linearly independent components (as representations of the linear group) in each part. 
Every part adopts a particular axisymmetric form, which can be deduced from the examination of the decomposition 
of the linearly independent components under rotations. 

This variety of components of the "viscosity tensor" is reflected in the various effective forces that arise from them. 
Some of these are longitudinal, that is, they are the gradient of a potential, and therefore do not contribute to the 
transverse rotating fluid equation. However, they arise from the turbulent state and contribute to the equation that 
determines the equilibrium state, so they may have practical relevance. Already at first order in f2, we have the 
potential fi • uj, which reminds us of the spin-orbit coupling of atomic physics. At higher orders in fi, we find more 
complicated potentials. 

Although the most general four-rank "viscosity tensor" includes terms that lead to the stress coupling to vorticity 
and not conserving angular momentum, one may wonder if they are really necessary. If we take the criterium of having 
the most general axisymmetric transverse force, we could remove redundant coefficients in Eq. (18): this equation 
has nine terms but includes the 14 coefficients of the traceless viscosity tensor. For example, we could remove all the 
coefficients belonging to (, and the couple C3, c 4 (belonging to x) or the couple c' 3 , c' 4 (belonging to £), but not c 2 or c 2 . 
Hence, we conclude that the only part of the viscosity tensor that can be neglected is (, but \ and £ must be present. 
So we still have that the stress couples to vorticity and does not conserve angular momentum (in the rotating frame). 

Finally, in our analysis of the four-rank tensors defined in some approaches to non-rotating anisotropic turbulence, we 
have seen there are similarities with our viscosity tensor in their definition and, therefore, that the respective definitions 
can be connected. However, the lack of specification of a quantity that breaks rotation invariance precludes that the 
actual values of the tensors corresponding to non-rotating turbulence can be anisotropic. It is possible, nevertheless, to 
provide such a symmetry-breaking quantity: for example, an anisotropic noise spectral correlation D mn . In particular, 
one can introduce an axisymmetric D mn by postulating the presence of a global vector (of unknown origin) , as in the 
perturbative approach of Ref. [9]. If this symmetry-breaking vector were axial instead of polar, the four-rank tensor 
Vijmn of section V should have the same form that our T]ij mn . 
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RESOLUTION OF THE GENERAL FOUR-RANK TENSOR BY THE SYMMETRY OF PAIRS OF 

INDICES 

Let us work out first the resolution of the general four-rank tensor T^ mn into a sum of tensors of definite symmetry 
type given by standard Young tableaux. Young tableaux indicate certain symmetry operations performed on the 
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indices [14]. We can consider the general four-rank tensor as a tensorial product of four vectors and, therefore, write 
its resolution as the Clebsch-Gordan decomposition for the linear group GL(3) of the corresponding direct product: 
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(40) 



The dimensions of the GL(3) representations on the right-hand side are: 15 for the totally symmetric representation, 
15 for the next mixed symmetry representation, 6 for the following mixed symmetry representation, and 3 for the last 
mixed symmetry representation. Therefore, wc have 81 = 15 + 3x15 + 2x6 + 3x3. 

We intend to show the correspondence of the preceding tensorial representations with the tensorial representa- 
tions selected according to the symmetry relative to pairs of indices. These are constructed as direct product of 
representations: 



U_LU 



\J ® — , ^ (§) — 



(41) 



The corresponding dimensions are: 36 for the symmetric-symmetric, 18 for both the antisymmetric-symmetric and 
the symmetric-antisymmetric, and 9 for the antisymmetric-antisymmetric. The symmetric-symmetric tensor and the 
antisymmetric- antisymmetric tensor can be further resolved into pair-symmetric and pair-antisymmetric components. 
Moreover, most of the six resulting representatins are still reducible. To find the irreducible representations, we will 
take advantage of the Clebsch-Gordan decomposition given by Eq. (40), superimposing on it the symmetry relative 
to pairs of indices. 

To have symmetry in the first pair, say, we may just symmetrize the general tensor (40) over indices ij. This 
operation inmediatly removes the fourth, sixth, eighth and ninth representations, which are antisymmetric in ij. In 
other words, given that 



\i\3\ 



(42) 



we can resolve the general tensor (40) into ij-symmetric and ij-antisymmetric parts. The former is 
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(43) 



where the right-hand side of the latter equation is the result of symmetrization over indices ij (unnecessary in the 
totally symmetric representation). The r/-antisymmetric part is more complicated but we will not need it. 
We can further resolve Eq. (43) by symmetrizing or antisymmctrizing over the remaining pair of indices: 
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Let us analyze Eq. (44). Using the definitions of the Young tableaux, we compute 

Tijmn ~i~ Tjimn ~i~ ^ijnm ~t~ ^jinm ^mnij ^mnji ^nmij ^nrnji ■> (^^) 

(ij),(mn) 
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rri rri rri rri rri rri rri rri 
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(47) 



If we denote the components of the tensor with symmetry by pairs [given by the left-hand side of Eq. (44)] as 

Sijmn = Tijran "t - Tji mn -\~ Tij nrn -\- Tji nm , (48) 

we can write Eq. (46) as 



and Eq. (47) as 



15 ^ijmn ~ Sijmn S mn ij 



G^ijmn — ^ijmn ~t~ S mn ij ^ [^imjn ~t~ Si nm j + S m ji n + <SVijmi] , 



where the left subscript indicates the dimension of the irreducible linear representations. On the other hand, the 
totally symmetric representation reads 



IhSijmn Sij mn + S mn ij + Sj m i n + Si m j n + Si n j m + >Sj r 

Now, let us analyze Eq. (45). We compute 
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(51) 
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If we denote the components of the tensor with symmetry in the first pair and antisymmetry in the second [given by 
the left-hand side of Eq. (45)] as 



we can write Eq. (50) as 



and Eq. (51) as 



SAij mn — Tij mn -\- Tji mn Tij nrn Tji nrai (52) 



lsSAijmn — SAij mn -\- ^ [SAimjn + SAi n7n j + SA m ji n + SAjijmi] 



In analogy with Eq. (52), we define 

ASij mn Tij mn Tji mn -\- Tij nm Tji nm , (53) 

and we have analogous definitions for the irreducible linear representations 15 ASij mn and 3ASij mn . 
We must also define the tensor with antisymmetry by pairs 

Aijmn Tij rnn Tji mn Tij nm -\- Tji nrn , (54) 

which can be divided into pair-symmetric and pair-antisymmetric tensors as 

6Aij mn Aij mn + A mn ij , (55) 
zAijvun Aij mn A mn ij . (56) 

Note that qA exactly corrsponds to the second square Young tableau of Eq. (40). 

Of course, many of the preceding irreducible representations of the group GL(3) are reducible with respect to 
its rotation subgroup 0(3), since the metric 5ij allows one to extract lower-rank representations by contracting 
indices. For example, the totally symmetric tensor contains the representations J — 4, 2,0 [14]. The six-dimensional 
representation given by qS or qA contain J = 2, 0. The 15-dimensional representations given by 15S" or 15SA contain 
J = 3, 2, 1. The irreducible representations of the group 0(3) are the symmetric traceless tensors [14]. Therefore, it 
must be possible to express each of the previous representations in terms of four-rank tensors as symmetric traceless 
tensors (of lower rank). For example, for J = 3, from i 5 S" we obtain T tjn = eu m S- jmn + ej im S' ilmn + e nim S' ijml ; 
for J = 3, from 15SA we obtain T p ij — e pmn SAij mn + timn SAj pmn + ej mn SA p i mn . Note that, in both cases, the 
components of the J = 2 representations vanish and do not contribute to the third-rank tensor, while the remaining 
J = 1 representation must be removed from this tensor by imposing that it be traceless. 
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